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The space mission CoRoT provides us with a large amount of high-duty cycle long-duration observations. Mode fitting 
has proven to be efficient for the complete and detailed analysis of the oscillation pattern, but remains time consuming. 
Furthermore, the photometric background due to granulation severely complicates the analysis. Therefore, we attempt 
to provide an alternative to mode fitting, for the determination of large separations. With the envelope autocorrelation 
function and a dedicated filter, it is possible to measure the variation of the large separation independently for the ridges 
with even and odd degrees. The method appears to be as accurate as the mode fitting. It can be very easily implemented 
and is very rapid. 
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1 Introduction 

A new era in asteroseismology has been opened by the mis- 
sion CoRoT, that provides us with high-quality photomet- 
ric time series. The Fourier spectra have revealed a large 
number of solar-like oscillations in main-sequence stars and 
subgiants, plus oscillations in thousands of giants. In com- 
parison to previous ground-based observations, observation 
duration is long, with a duty cycle of the order of 90 %. The 



signal in the time series. For rapidity, the autocorrelation of 
the time series is performed as the Fourier spectrum of its 
Fourier spectrum. For efficiency and detailed analysis, a fil- 
ter helps to select a given frequency range in the Fourier 
spectrum. MA09 have analyzed the performance of the de- 



tection of the large separation with the EACF. Mosser & Appourchaux (2010) 
have presented the structure of the automated pipeline based 
on the EACF. 



The method has proven to be efficient at low signal-to- 
ground-based network dedicated to Procyon ('SentoftetaTlO^ise ratio (Mos ser et al.^009] ). A positive detection, with 



has reached a similar duty cycle, but only for the central 
10 days of the network observation. 

Both the analysis of these solar-like oscillations and the 
scientific output crucially depend on the quality of the data. 
Automated pipelines based on different methods have been 



> ■ setup jMathur et al. 2010llHekker et al. 2009l|Huber et al. 2009t egrees of the modes 



the maximum amplitude of the EACF above the threshold 
level, gives at least the measurement of the mean value of 
the large separation. At larger signal-to-noise ratio (SNR), 
the method provides the variation of the large separation. 
For high-quality data, it gives also a test for identifying the 



In order to extract precise eigenfrequencies, modes fitting is 
required, as made e.g. in |Benomar et al. (2009)| The perfor- 
mance of this time-consuming method appears to be lim- 
ited by the background level. Therefore, we investigate an 
alternative approach. We focus our analysis on the mea- 
surement of the large separation and propose to extend the 



method presented in jMosser & Appourchaux (2009) here- 
after MA09. 



Roxburgh & Vorontsov (2006) have presented how to an 



alyze solar-Uke oscillations with the autocorrelation enve- 
lope function (EACF). This function searches directly for 
the signature of the large separation in the autocoiTelation 
of the time series, without any prior on the oscillation ex- 
cess power The EACF measures directly the acoustic ra- 
dius of the star The delay of coiTelation being shorter than 
the mode lifetimes, the correlation is able to stack coherent 
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Deriving the variation of the large separation with fre- 
quency is highly interesting for stellar modeling (IDeheuvels et al. 2010bl l. 
Contrary to eigenfrequencies, the large separations are much 
less sensitive to the influence of the upper atmospheric lay- 
ers. However, as shown by many studies (e.g. lProvost et al. 19931 1. 
the large separation variation /S.v{v) may depend on the de- 
gree of the mode. A method for deriving lS.v{v) for indi- 
vidual degree is highly desirable. We will consider the mea- 
surement of Ai^cvcn and Ai/odd, corresponding to the even 
and odd values of the degree. [Bedding & Kjeldsen (2010)| 
have proposed to make use of the frequencies of the ridge 
centroids in order to define Ai/cvcn and Az^odd- In this pa- 
per, we propose a more direct approach, based on the EACF 
but with a dedicated filtering. 

The method is presented in Section|2] We present its co- 
herence: it first gives the mean value (Az^) of the large sep- 
aration, then the variation /S.v{v), then the values At'oven 
and Ai^odd obtained for the even and odd ridges separately. 
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Different results are presented and discussed in Sect. |3] We 
show that we easily reproduce results based on the time con- 
suming mode fitting. Section|4]is devoted to conclusions. 

2 Method 

We intend to perform an analysis of the even and odd ridges 
separately. Therefore, the basic concept of the method con- 
sists in introducing a comb structure in the filter of the EACF. 
We present the recipe used for extracting such an infor- 
mation without a priori. We use the EACF to first define 
the comb and then derive the measurement of Ai^cvcn and 
At'odd- We do not base the analysis on the solar case, since 
the Sun does not present significantly different Ai/cvcn and 
Ai/odd- Instead, we consider the case of HD 49385, a G- 
type subgiant observed with CoRoT that shows significantly 
different even and odd ridges dPeheuvels et al. 2010al l: the 
large separations between radial modes are about 56 fjHz, 
whereas for ^ = 1 modes they can be as low as 50 /LtHz. 

2.1 Analysis of the comb structure 

MA09 have determined that the EACF with a filter width 
equal to the full-width at half-maximum (FWHM) of the 
excess power envelope centered on the frequency i^max of 
maximum amplitude gives the mean value ( Ai^) of the large 
separation. With a narrower filter, with a typical FWHM of 
2(Azy), one can extract the variation of the large separation 
with frequency Ai'{i'). In order to analyse independently 
the even and odd ridges, we introduce a comb structure in 
the filter J-av- 




(mHz) 



Fig. 1 Aj/| 1 1 (i/) for HD 49385 based on 5 different large 
separations varying from 54 to 58 /iHz. The inset shows the 
filter J^Aiy for (Aj/) — 56 /iHz and Siyp = 2.5 (A;/). 
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Fig. 2 Same as Fig. [T] with the superimposition of the 
estimates of the eigenfrequencies and large separations ob- 
tained with Eq.[3] The filter is based on the mean large sep- 
aration (Ai/); its FWHM is changed compared to Fig.[T]in 
order to show its reduced influence. 



for \v'\ > Svp 



where v' is the frequency shift to a central frequency and 
Si^F is the FWHM of the Hanning filter We cafl Ai/|||(:/) 
the function Ah'{i') obtained from the autocorrelation with 
such a filter (Fig.[T]). 

The definition of Aj/| | | (t^) has the disadvantage to intro- 
duce the mean value (Aiy) of the large separation into the 
analysis. Subsequently, it may be suspected to bias the re- 
sult. We have shown that it is not the case, by comparing the 
values of Ai^i 1 1 (i^) obtained with different comb filters based 
on different large separations (Fig. [Til. We have superim- 
posed the results for 5 different values of the comb period: 
(Av), {Av) ± 1 /xHz and {Au) ± 2 //Hz. The ±2 ^Hz shift 
represents a relative variation of about 4 %, much greater 
than the precision of the measurement of {Av). Since the 
different values of Av\ \\(y), even when obtained with strongly 
biased filters, are in fact very similar, we conclude that tak- 
ing a filter based on the mean value of the large separation 
is adequate. 



2.2 Ridges 

The function Av\\\{v) previously defined is continuous. In 
order to derive Av^vcn and Ai^odd, we have to identify the 
even and odd ridges and to locate the eigenfrequencies inde- 
pendently of any mode fitting. The ridge identification was 



already performed in MA09 or by Bedding & Kjeldsen (2010) 
so that we are left with the estimate of the eigenfrequencies. 

In order to perform this, we simply use the EACF to 
measure the first order variation of the large separation: 

Au ~ {Au) +an' (2) 

with n' ~ n — riQ + £/2 and no the radial order associ- 
ated with the nearest radial eigenfrequency to t'max- From 
Eq. (|2]), we can derive proxies of the ^ = and £ = 1 eigen- 
frequencies: 



l^n,0 ^Vo+ Al^ovon n' + ttQ 



(3) 
(4) 
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Fig. 3 Ai^ovon and Avodd for HD 49385. The solid 
lines (respectively dashed lines) correspond to the even 
(reps, odd) ridge; l-cr error bars are given with the same 
linestyles. We have superimposed the values presented in 
IDeheuvels et al. 2010al derived from the mode fitting: dia- 
monds correspond to £ — 0, triangles to i' = 1. Error bars of 
the fitted values, not shown, are equivalent to the error bars 
derived from the EACF. 



All parameters are direct outputs of the EACF analysis, ex- 
cept the constant terms vq and i^i. Their precise values, 
close to I'niax, are derived from the minimization of the 
residuals between the observed and estimated eigenfrequen- 
cies. 

Figure|2]shows the values of the large separation Ai/| 1 1 (i^) 
estimated at the frequencies determined by Eq. (O. We note 
that Avcvcn and Ai'odd correspond to local extrema of the 
function Ai/|||(i^). This justifies a posteriori that a precise 
location of the eigenfrequencies is not necessary to derive 
Ai'ovcn and Ai^odd- We can then compare Ai^cvcn and At'odd 
to their fitted values, namely the values derived by the mode 
fitting. 



2.3 Performance 



The performance of the method, expressed by the error bars 
on Ai^cvcn and Az^odd, depends basically on the scaling of 
the noise (Eqs. 3-6 of MA09). Compared to the results pre- 
sented in MA09, the error bars are essentially multiplied by 
a factor 2 since the comb pattern divides by 2 the efficient 
number of points in the filter. 

Tests with different filters have shown that rectangular 
or triangular shapes induce too many discontinuities in the 
function Ai^m (i^). We therefore prefer to consider Hanning 
filters, as in MA09, with a FWHM varying from 2 to 4 
(Ai^). Smaller widths give the finest resolution, but require 
a high signal-to-noise ratio, indicated by high values of the 
EACF 




Fig. 4 Same figure as Fig. [3] comparing Az/ovcn and 
Ai'odd for HD 49933 with fitted values determined by 
IBenomar et al. (2009)| 

HD181420 




iHz) 

Fig. 5 Same figure as Fig. [3] comparing Az/ovcn and 
Avodd for HD 181420 with fitted values determined by 
[Barban et al. (2009)| Data are more noisy, but the agreement 
is as clear as for oscillation spectra with higher signal-to- 
noise ratio. 

Procyon 




I' (mHz) 

Fig. 6 Same figure as Fig. [3] comparing Az/cvcn and 
Avodd for Procyon with the centroids of the ridges (Fig. 



1 la of Bedding et al. 2010| . A larger filter has been used to 
compare our results to the smoothed determination provided 
by these centroids. 
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Fig. 7 Same figure as Fig. [3] comparing Ai^cvcn and 
Afc'odd for HD 181906 with fitted values determined by 
IGarcia et al. 20091 
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Fig. 8 Same figure as Fig. [3] comparing Az/ovon and 
Aj^odd for the star KIC 01 1026764 observed with Kepler, 



3 Discussion 

3.1 A few results 

We have compared for different stars Ai^ovcn and At^odd 
obtained with this method or derived from the mode fitting. 
The agreement lies within the error bars (Fig.[3]to|6l). We 
note that the error bars obtained with the EACF are equiva- 
lent to the l-cr error bars of the fitted eigenfrequencies. For 
clarity, we do not reproduce them in the figure. 

For HD 49385 (Fig. O, the EACF analysis is clearly 
able to disentangle the very different values of Ai/even and 
Ai^odd due to avoided crossings (Deheuvels et al. 2010al) . 
For HD 49933, the method exhibits both the significant mod- 
ulation of Aiy{h') and the difference between Ai/cvcn and 
At^odd- 

3.2 At low signal-to-ratio 

For HD 1 8 1420, observed with a lower signal-to-noise ratio, 
the agreement remains remarkable (Fig. |5]l. The values of 
Ai/cvcn and Ai/odd extracted from the EACF have smoother 
variation than given by the mode fitting. Due to the value 
of error bars, one cannot decide if the difference is due to 
the smoothing provided by the filter of the EACF or to spu- 
rious variation in the fitted values. On the other hand, the 
EACF provides values for Afcvon and Ai'odd in a larger fre- 
quency range. By the way, we confirm the identification of 
the degree of the ridge proposed by MA09, corresponding 
to scenario 1 of |Barban et al. (2009)| 

It is possible to use a larger filter in order to analyze 
low SNR light curves. We have performed the analysis of 
Avi 1 1 (i/) for the CoRoT target HD 1 8 1906 (IGarcia et al. 20091) . 
The difference between Ai'cvcn and Ai^odd (Fig-Hl is clearly 
shown, but the method does not provide a way to identify 
the ridges. We believe that the clean measurement of Aj/ovcn 
and Ai^odd may be useful for comparison with a forthcom- 
ing modeling of this star 



with fitted values determined by Chaplin et al. (2010)[ The 
function Ai/|||(i') has been superimposed (dotted lines), 
with discontinuities in the frequency range where the reg- 
ularity of the oscillation pattern of ^ = 1 modes is broken. 



HD49933 




Fig. 9 A2i'cvon and A2i'odd for HD 49933, presented as 
Av in Fig.Oand compared to the fitted second differences 
derived from |Benomar et al. (2009)| 



HD52265 




(mHz) 



Fig. 10 Same as Fig. |9] for HD 52265, with comparison 
to the results from.Ballot et al. 2010t 
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3.3 Avoided crossings 

The presence of avoided crossings may complicate the mea- 
surement of Az^i 1 1 (i^), since they perturb the regular Tassoul- 
like pattern. We have tested this effect on the evolved star 



KIC 1 1026764 observed by Kepler and analyzed by ChapUn et 



This method provides an alternative to mode fitting. It 
does not give the precise eigenfrequencies, amplitudes and 
lifetimes of a solar-like oscillation pattern. However, with 
the variation Ai^(i^) of the large separation, with Afovon and 
Aii:ftdd.,-aiid with accurate proxies of the eigenfrequencies. 



This star shows mixed modes, with a complex pattern in the 
frequency range [0.88, 0.95 mHz]. We show that the method 
is sensitive to the irregularity of the oscillation pattern due to 
avoided crossings. The discontinuities in /S.v\\\{v) indicate 
clearly the presence of the mixed modes (Fig. [8]). Instead 
of pointing the perturbed value Aj^avoid of the large sepa- 
ration at the avoided-crossing frequency, the method gives 

2Az>'cvcn Az^avoid- 

3.4 Second difference 

Measuring the large separations gives access to further im- 
portant asteroseismic typical frequencies. The so-called sec- 
ond differences 

^2VnJ = Vn+l4~2Vn,l + Vn-l,l — Al/„+i^^ — Az^n^f (5) 

can be estimated directly. These second differences are of 
great help to investigate structure discontinuities, as for ex- 
ample the depression in Fi observed in the second helium 
ionization zone, that can be used to determine the enve- 
lope helium abundance of low-mass main-sequence stars 
dBasu et al. 2004i 

Since the method makes it possible to measure large 
separations depending on the parity of the degree, we can 
obtain odd and even values of the second differences, re- 
spectively defined as: 

A2t'ovcn/odd = Al'evcn/odd('^ + 1) ~ Al'ovcn/odd('^)(6) 

From the results presented in Fig. 21 we derive these even 
and odd second differences for HD 49933 (Fig.|9]l. They show 
a very close agreement with the values derived from the fit- 
ted individual eigenfrequencies. Fig. [TO] presents the varia- 
tion of /S.2V for HD 52265, another star observed by CoRoT 
jBallot et al. 201 Oi l, with again a close agreement. In both 
cases, we remark a phase shift between the modulation of 
A2i'cvcn and /S.2Vodd- This reinforces the interest to mea- 
sure the even and odd second differences according to Eq. 

m 



cess to a large part of the asteroseismic investiga- 
tion. Its advantage, compared to mode fitting, consists in its 
rapidity. 
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4 Conclusion 

We have shown that the envelope autocorrelation function 



presented by Mosser & Appourchaux (2009) and used with 



an ad-hoc filter is able to give the measure of the large 
separation of the even and odd ridges. The method is self- 
consistent: the EACF first provides the mean value of the 
large separation, then the function Av{v) derived with a 
narrower filter, and finally the functions Az/even and Ai^odd 
derived with a comb filter. From these values, we can derive 
the second differences. Error bars provided by the EACF are 
of the same order as error bars of the mode fitting. 
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